Fuel cell bipolar plates with flow uniformity

ABSTRACT

A bipolar plate with an enhanced fluid flow field design is provided for improved flow uniformity. The bipolar plate includes an inlet, an outlet, and a flow field having a pattern defining a plurality of microchannels configured to provide fluid communication between the inlet and the outlet. The pattern includes a plurality of channels with discrete areas of discontinuity to direct fluid from the inlet to the outlet. The pattern is designed using an inverse permeability field and is based on a reaction-diffusion algorithm to model channel spacing, thereby providing a variable pitch microchannel pattern. In various aspects, Gray-Scott reaction-diffusion equations may be used to obtain an anisotropic microchannel layout. Methods may include optimizing a porous media model domain for one of: a minimum flow resistance across the domain, a uniform velocity across a direction of the domain, and a minimum velocity across a direction of the domain.

TECHNICAL FIELD

The present disclosure generally relates to fuel cell bipolar plates and, more particularly, to designs for tailoring and customizing the flow fields of bipolar plates to control a fluid flow pattern and flow uniformity.

BACKGROUND

The background description provided herein is for the purpose of generally presenting the context of the disclosure. Work of the presently named inventors, to the extent it may be described in this background section, as well as aspects of the description that may not otherwise qualify as prior art at the time of filing, are neither expressly nor impliedly admitted as prior art against the present technology.

Bipolar plates useful with fuel cell technology are designed to permit the transfer of air or fuel from an inlet, through a flow field, and to an outlet. The flow fields of the bipolar plates may include channels commonly formed by machining or stamping techniques. The efficiency of the chemical reaction within the fuel cell is based, in part, on the proper flow and distribution of air and fuel fluid flow streams. The fluid distribution is, in turn, based on the channel network and design. A uniform fluid flow is important to both performance and cost. Various designs of flow fields for bipolar plates often include portions with straight channels, and portions with channel designs that are manually optimized by using computational fluid dynamics in order to determine which portions of a flow field should be solid, and which portions should be provided for the fluid transfer. The manual optimization is both time consuming and costly.

Accordingly, it would be desirable to provide improved, cost effective bipolar plate designs and methods of manufacture that include an optimized fluid flow field distribution that can be customized.

SUMMARY

This section provides a general summary of the disclosure and is not a comprehensive disclosure of its full scope or all of its features.

In various aspects, the present teachings provide a method for designing a microchannel layout for a flow field of a bipolar plate. The method may include defining a fluid flow optimization domain with boundary conditions and loads. The fluid flow optimization domain has an x-axis defining a longitudinal direction from an inlet leading to an outlet window of the flow field, a y-axis defining a transverse direction, and a z-axis perpendicular to the x-axis and the y-axis. The method may include using a gradient-based algorithm together with computational fluid dynamics to optimize a porous media model domain. This may include setting a minimum inverse permeability to a non-zero value and obtaining a grayscale design and fluid velocity field. The method may then use Gray-Scott reaction diffusion equations with the grayscale design and fluid velocity field in order to obtain a microchannel layout with a plurality of channels. The method includes providing the plurality of channels with discrete areas of discontinuity and incorporating the microchannel layout as a pattern for an inlet region of the flow field of the bipolar plate. In various aspects, the discrete areas of discontinuity are disposed along spaced apart contour lines based on a predetermined performance parameter of the bipolar plate

In other aspects, the present teachings provide a method for designing an anisotropic microchannel layout for a flow field of a bipolar plate. The method may include defining a fluid flow optimization domain with boundary conditions and loads. The method may include using a gradient-based algorithm together with computational fluid dynamics to optimize a porous media model domain. This may include setting a minimum inverse permeability to a non-zero value and obtaining a grayscale design and fluid velocity field. The method may then use the Gray-Scott reaction diffusion equations and an anisotropic diffusion tensor along with the grayscale design and fluid velocity field and to generate an anisotropic microchannel layout having a connected line and space pattern including a plurality of channels. The plurality of channels are provided with discrete areas of discontinuity. The anisotropic microchannel layout is incorporated as a pattern for an inlet region of the flow field of the bipolar plate to provide flow uniformity.

In still other aspects, the present teachings provide a bipolar plate for a fuel cell including an inlet, an outlet, and a flow field. The flow field includes an inlet region having a pattern defining a plurality of anisotropic microchannels configured to provide fluid communication between the inlet and the outlet. The anisotropic microchannels are defined having discrete areas of discontinuity. The pattern of the flow field is designed using an inverse permeability field and solving Gray-Scott reaction-diffusion equations with an anisotropic diffusion tensor to obtain an anisotropic microchannel layout with channel spacing based on effective medium theory. In various aspects, the discrete areas of discontinuity are disposed along spaced apart contour lines based on a predetermined performance parameter of the bipolar plate selected from at least one of pressure flow and velocity flow.

Further areas of applicability and various methods of enhancing the above technology will become apparent from the description provided herein. The description and specific examples in this summary are intended for purposes of illustration only and are not intended to limit the scope of the present disclosure.

BRIEF DESCRIPTION OF THE DRAWINGS

The patent or application file contains at least one drawing executed in color. Copies of this patent or patent application publication with color drawings will be provided by the Office upon request and payment of the necessary fee.

The present teachings will become more fully understood from the detailed description and the accompanying drawings, wherein:

FIG. 1 illustrates a schematic plan view of an exemplary bipolar plate assembly for a fuel cell useful according to various aspects of the present technology;

FIG. 2 is a magnified view of a portion of the bipolar plate assembly of FIG. 1 represented by the shape labeled with reference number 2;

FIG. 3 is an alternate view of the portion of the bipolar plate assembly of FIG. 2 provided with a simplified inlet design and illustrating a magnified portion of a flow field outlet window region;

FIG. 4 illustrates an exemplary schematic flow chart for designing a pattern for an inlet region of a flow field;

FIG. 5 is a flowchart of an exemplary gradient based algorithm using method of moving asymptotes (MMA) or globally convergent method of moving asymptotes (GCMMA) looped in with the fluid flow CFD to optimize the domain;

FIG. 6 provides a schematic illustration a fictitious design variable of density γ that may be introduced to relax the topology optimization problem and provide a numerical framework to establish flow optimization;

FIG. 7 provides an illustration of a flow field with the inverse permeability κ_(min) being equal to zero, and represents a strictly black & white (0, solid and 1, fluid) design;

FIG. 8 provides an illustration of an exemplary grayscale design representing a porous medium with the inverse permeability κ_(min) being equal to one third of the maximum, for example, κ_(min)=0.333 κ_(max);

FIG. 9 illustrates an exemplary method flow chart for designing a pattern for an inlet region of a flow field using an anisotropic diffusion tensor;

FIG. 10 illustrates an exemplary schematic generally summarizing the methods for designing a pattern for an inlet region of a flow field;

FIG. 11 illustrates an exemplary pressure field obtained from an optimized porous media inlet flow field model;

FIG. 12 illustrates a line graph of a velocity profile associated with the pressure field of the optimized porous media inlet flow field model of FIG. 11;

FIG. 13 illustrates a resulting pressure field obtained from a continuous channel design of the optimized porous media inlet flow field model of FIG. 11;

FIG. 14 illustrates a line graph of a velocity profile associated with the pressure field of the continuous channel inlet flow field design of FIG. 13;

FIG. 15 illustrates a resulting pressure field obtained from a discontinuous channel design of the optimized porous media inlet flow field model of FIG. 11 having a substantially uniform inlet feed rate;

FIG. 16 illustrates an exemplary design of a flow field with a discontinuous channel inlet flow field design according to another aspect of the present technology and having a non-uniform inlet feed rate;

FIG. 17 illustrates an exemplary velocity field obtained from an optimized porous media inlet flow field model;

FIG. 18 illustrates a resulting velocity field obtained from a discontinuous channel design model of the optimized porous media inlet flow field model of FIG. 17;

FIG. 19 is a graph that shows the results of a comparison between a discontinuous channel inlet velocity field and a target specification velocity field; and

FIG. 20 illustrates a cross-sectional view of an exemplary set of channels and walls that can be made according to various manufacturing techniques.

It should be noted that the figures set forth herein are intended to exemplify the general characteristics of the methods, algorithms, and devices among those of the present technology, for the purpose of the description of certain aspects. These figures may not precisely reflect the characteristics of any given aspect and are not necessarily intended to define or limit specific embodiments within the scope of this technology. Further, certain aspects may incorporate features from a combination of figures.

DETAILED DESCRIPTION

The present technology generally teaches an optimized bipolar plate structure and/or assembly for a fuel cell, along with methods of design and manufacturing. The bipolar plates used in the structure and/or assembly can be provided with a flow field that can be designed with a plurality of channels, also referred to as microchannels, to provide the bipolar plate with a tailored and customized design. The channels of the present technology are defined having discrete areas of discontinuity that provide an increased uniformity in fluid flow. The pattern of the flow field may be designed using an inverse permeability field and solving Gray-Scott reaction-diffusion equations with an anisotropic diffusion tensor to obtain an anisotropic microchannel layout with channel spacing based on effective medium theory. In various aspects, at least a portion of the discrete areas of discontinuity are disposed along spaced apart contour lines that may be based on a predetermined performance parameter of the bipolar plate. For example, the predetermined performance parameter may be at least one of pressure flow and velocity flow.

In various aspects, methods for designing a microchannel layout for a flow field of a bipolar plate include defining a fluid flow optimization domain with boundary conditions and loads. The methods may include performing an interpolation process to design a flow distribution network. Using a gradient-based topology optimization algorithm, looped in together with fluid flow computational fluid dynamics, the domain is then optimized for minimum flow resistance (with optional mass flow rate constraints at fluid outlet boundaries) using a modified Darcy flow equation for porous media. In other aspects, the domain optimization can be based on obtaining a uniform velocity across a direction of the domain, or a minimum velocity across a direction of the domain. In contrast to various topology optimization techniques, a non 0-1 (solid-fluid) design is obtained with grayscale throughout the domain. For example, a material interpolation approach is used where the method includes setting the minimum inverse permeability to a non-zero value (with the maximum permeability representing a porous medium), and obtaining a grayscale design and fluid velocity field. Using Gray-Scott reaction diffusion equations with the grayscale design and fluid velocity field, the method includes obtaining a microchannel layout. The channel spacing can be based on effective medium theory for flow through porous media. In various aspects, the permeability of the optimized design space is linked to the microchannel design layout through an anisotropic definition of the diffusion coefficients of the Gray-Scott reaction diffusion equations in order to obtain an anisotropic microchannel layout. The use of present technology eliminates the explicit modeling of the microchannels. Instead, it derives a channel design using a computationally efficient post-processing technique based on the Gray-Scott equations and the velocity flow field through the porous medium. The discontinuous channel design is able to recover the porous model pressure field distribution with a high degree of accuracy, and also provides a match between models for an optimized velocity field. For example, in various aspects, the discontinuous channel design is able to meet target velocity flow rate specifications with less than about 5% of an average outlet flow rate variation.

The optimized design for the microchannel layout is then incorporated as a pattern for the flow field of the bipolar plate. For example, fabrication techniques can be used, in combination with the pattern, to create a network of functionally graded microchannels, a porous media, or a hybrid combination thereof. As will be described in more detail below, additive manufacturing techniques can be used where the flow field, or portion thereof, is fabricated using laser sintering, e-beam melting, and/or binder jet technology that may use a combination of printing and curing. Electroplating techniques may also be used for fabricating a functionally graded porous medium. A commercially available permeable material or functionally graded porous medium can be used. In other aspects, the porous medium can be custom made. In some examples, metal inverse opal technology can be used with a suitable metal, such as copper, nickel, or titanium. The functionally graded porous medium is generally obtained by templating polystyrene particles of different sizes, and/or controlling a sintering time for different templated regions. Conventional metal machining and/or stamping techniques can be used in combination with the above techniques to ultimately provide the bipolar plates.

FIG. 1 illustrates a schematic plan view of an exemplary bipolar plate assembly 20 that may be used with a fuel cell (not shown), such as a PEM fuel cell, according to various aspects of the present technology. Bipolar plate assemblies 20 are important components of PEM fuel cells, and have commonly been made of graphite or metals, such as stainless steel, aluminum, and titanium; composite bipolar plates have also been used. Bipolar plate assemblies 20 typically supply fuel and oxidant to reactive sites adjacent the membrane, remove reaction products, collect current produced by the internal reaction, and provide the necessary mechanical support for a plurality of fuel cells in a fuel cell stack.

By way of background, bipolar plate assemblies may include an anode flow field plate and a cathode flow field plate that have been bonded or otherwise appropriately sealed together as an assembly. In certain instances, they may be bonded to form a sealed coolant flow field between the plates, as commonly employed in the art. Various transition channels, ports, ducts, and other features involving all three operating fluids (i.e. fuel, oxidant, and coolant) may also appear on the inactive side and other inactive areas of bipolar plate assemblies. The operating fluids may be provided under significant pressure, thus it is important that all of the features in the plates are appropriately sealed to prevent leaks between the fluids and to the external environment. Another feature for bipolar plate assemblies is that there is a satisfactory electrical connection between the two plates because a substantial current generated by the fuel cell stack must be able to pass between the two plates.

The plates making up the assembly may optionally be metallic and are typically produced by stamping the desired features into sheets of appropriate metal materials (e.g. certain corrosion resistant stainless steels). Two or more stamped sheets are then typically welded or clamped together with gaskets so as to appropriately seal all the fluid passages from each other and from the external environment. Additional welds may be provided to enhance the ability of the assembly to carry electrical current, particularly opposite the active areas of the plates. Metallic plates may however be bonded and sealed together using adhesives or the aforementioned gaskets. Corrosion resistant coatings are also often applied before or after assembly.

In other aspects, the present technology provides monolithic and/or hybrid structure bipolar plate assemblies. As used herein, the term “monolithic” means a single, unitary component that is intractably indivisible once formed. While there may be an interface region where two or more components are joined together by bonding or some type of fusion, once the materials are bonded/fused to one another, the materials do not typically separate from one another as may occur in other layered substrates. The bipolar plates may include porous materials and functionally graded materials as a base structure, and/or may include portions that include a material created by additive manufacturing techniques, electroplating techniques, and the like. As used herein, the term “hybrid” structure means a structure that includes more than one type of material or structure, for example, one non-limiting aspect may include a first portion including a porous medium, and a second portion with explicit microchannels, which may be created by additive manufacturing techniques, electroplating techniques, and the like. In other non-limiting aspects, a hybrid porous structured material may include a first portion/region that is non-porous, and a second portion/region that is porous with a certain permeability. The design and manufacturing of monolithic and hybrid structure types of bipolar plates and/or bipolar plate assemblies made with additive manufacturing techniques and electroplating techniques will be described in more detail below.

With renewed reference to FIG. 1, the structure of the bipolar plate assembly 20 generally includes at least one basic support substrate or bipolar plate structure 22 defining one or more inlet 24 and outlet 26, and a flow field 28 that cooperates with the bipolar plate assembly 20 providing fluid communication between a respective inlet 24 and outlet 26. The flow field 28 may be provided with one or more different regions having a respective shape or pattern that defines a plurality of channels 30, or microchannels, sized and configured to provide the appropriate fluid communication. While one or more of the regions may be provided with generally straight, parallel, and/or serpentine style flow channels, in various aspects, certain of the channels and/or portions of the flow field, such as the inlet region 32 and/or outlet region 34, may be shaped and oriented in order to specifically minimize flow resistance and to provide a minimal fluid power/pressure drop or to provide a uniform flow field across the plurality of microchannels. In addition to providing fluid communication between an inlet 24 and a respective outlet 26, the channels 30 may serve as a support structure for other components and the membrane and provide fluid transfer of the fuel and air to a respective anode diffusion layer and a cathode diffusion layer, which are then directed to the catalyst layers and reaction sites. As will be described in more detail below, various channels or portions thereof may be provided with discrete areas of discontinuity in the channel walls. This may provide for “branching” of the channels, which may operate similar to a manifold, allowing fluid flow through different branches, into and through several different openings.

FIG. 2 is a magnified view of a portion of the bipolar plate assembly 20 of FIG. 1 represented by the shape labeled with reference number 2. The portion shown in FIG. 2 includes an inlet 24 and the inlet region 32 of the flow field 28. As shown in FIG. 3, for ease of design, modeling, and explanation herein, the inlet 24 is simplified as having a circular shape that is then used for obtaining the domain and data for optimization described herein. The present technology generally focuses on the design of a flow pattern and microchannel designs of the inlet region 32 of the flow field 28, as shown in FIG. 2, so that it can be designed and configured with a plurality of microchannels and discrete areas of discontinuity to provide the inlet region 32 of the bipolar plate with a tailored and customized porosity and permeability that includes an improved flow uniformity. In broad terms, porosity is generally a measure of how much of a substrate is open space. This open space can be, for example, between grains or within cracks, cavities, and/or channels provided within or on a surface of the substrate. Permeability is generally a measure of the ease with which a fluid can move through the porous substrate.

FIG. 4 provides an exemplary schematic flow chart for designing a pattern for an inlet region of a flow field. The method may begin with calculating minimum and maximum permeability values from design constraints in order to obtain an effective permeability of a channel array. The fluid flow optimization domain, boundary conditions, and loads may be defined. Minimum and maximum values to bound the permeability field value may be used to execute topology optimization. As is known in the art, topology optimization is a mathematical method that optimizes a material layout within a specific design space, based on a particular set of loads, boundary conditions, and constraints, in order to maximize the performance of the specific system. In various aspects, a gradient-based algorithm, together with computational fluid dynamics (CFD), may be used to optimize the domain for minimum flow resistance. FIG. 5 is a flowchart of an exemplary gradient-based algorithm using method of moving asymptotes (MMA) or globally convergent method of moving asymptotes (GCMMA) looped in with the fluid flow CFD to optimize the domain. In this instance, the optimization determines where in a flow field design domain to place the channels and where to place the fluid. In various aspects, results obtained from certain topology optimizations can be directly manufactured using various additive manufacturing techniques. For the topology optimization, one general approach is to represent a structural configuration by a material presence at each point. As such, the present technology provides a flexible method for design of fluid flow manifolds under laminar flow considerations. With reference to FIG. 6, a fictitious design variable of density γ may be introduced to relax the problem and provide a numerical framework to establish flow optimization. FIG. 6 shows the various possibilities, with γ=1 representing a fluid state (essentially white in FIG. 6), γ=0 representing a solid or low-permeability solid, or solid (essentially black in FIG. 6). The grayscale coloring in FIG. 6 represents an intermediate state with 0<γ<1, somewhere between a fluid state and a solid state. Topology optimization is then applied to the solid-fluid manifold design problem.

Based on the incompressibility condition of ∇·u=0, where u is the fluid velocity, one can arrive at the Navier-Stokes equation for flow through a porous medium, as follows: ρ(u·∇u)=−∇P+∇{η[∇u+(∇u)^(T)]}−κ_(e)(γ)u where κ_(e) represents the effective inverse permeability. As a non-limiting example, the effective inverse permeability (convex) interpolation may then be provided as follows, with permeability as a function of γ, the design variable:

${\kappa_{e}(\gamma)} = {\kappa_{\min} + {\left( {\kappa_{\max} - \kappa_{\min}} \right)\frac{q\left( {1 - \gamma} \right)}{q + \gamma}}}$

Methods of designing a microchannel layout for the inlet region 32 of a flow field begin with defining a fluid flow optimization domain with boundary conditions and loads. A gradient-based algorithm may then be used, together with computational fluid dynamics (CFD), to optimize the domain for minimum flow resistance. FIG. 5 illustrates the use of finite element analysis to update design variables using an exemplary method of moving asymptotes (MMA) or a globally convergent method of moving asymptotes (GCMMA) algorithm that is looped in with the fluid flow CFD. In various aspects, the methods for minimizing the flow resistance include minimizing the fluid power and/or pressure drop. For the fluid inlet boundary conditions, a first option includes a specified velocity and a second option includes a zero pressure; other conditions such as a specified mass flow rate are possible, as well. For the fluid outlet boundary conditions, the first option includes a zero pressure, and the second option includes a specified velocity. Subject to a solid material volume constraint, for example 50%, the objective function F to minimize is based on the following relationship:

$F_{o} = {\int\limits_{\Omega}{\left\lbrack {{\frac{1}{2}\eta{\sum\limits_{i,j}\left( {\frac{\partial u_{i}}{\partial x_{j}} + \frac{\partial u_{j}}{\partial x_{i}}} \right)^{2}}} + {\sum\limits_{i}{{\kappa_{e}(\gamma)}u_{i}^{2}}}} \right\rbrack d\;\Omega}}$ The function generally represents the quantity that is being minimized for best performance, here that is the minimization of flow resistance across the structure for better efficiency. Minimization of flow resistance is important because it is directly related to pumping power. In various aspects, the design space Ω may be defined to provide regions that cannot be modified during the course of the optimization. Evaluation of the objective function can be performed solving differential equations using a finite element method or similar computational approach.

In other aspects, the methods may include using the gradient-based algorithm together with computational fluid dynamics to optimize the porous media model domain for a uniform flow velocity or a minimum flow velocity in a certain direction, for example, along an axis. With renewed reference to FIG. 1, the fluid flow optimization domain may include an x-axis defining a longitudinal direction from an inlet 24 leading to an outlet 26 of the flow field, a y-axis defining a transverse direction with respect to the x-axis, and a z-axis perpendicular to the both x-axis and the y-axis. In certain aspects, a flow velocity can be optimized for uniformity by taking into account a difference in the flow velocity in a direction of the x-axis is minimized according to a relationship of:

$\left\lbrack \frac{u - u_{{out}_{avg}}}{u_{{out}_{avg}}} \right\rbrack^{2}$ wherein u represents a local velocity in the direction of the x-axis exiting a channel of the outlet window, and u_(out avg) represents an optimized target velocity in the direction of the x-axis exiting a channel of the outlet window. In certain other aspects, the methods may include using the gradient-based algorithm together with computational fluid dynamics to optimize the porous media model domain for a minimum flow velocity in a direction of the y-axis. In this regard, the methods may include determining a difference in flow velocity in a direction of the y-axis that is minimized according to a relationship of: v², wherein v represents a local velocity in the direction of the y-axis.

FIGS. 7 and 8 illustrate the comparison between traditional topology optimization and porous media design. FIG. 7 provides the situation with the inverse permeability κ_(min) being equal to zero, and represents a strictly black & white (0, solid and 1, fluid) design, illustrating pure fluid flow from the inlet 24 shown with fluid streamlines, generally represented by reference number 36 that, as shown, exits via roughly 146 outlets (each outlet is not separately shown). One issue that arises with the configuration of FIG. 7, however, is that when there is a large open volume for a flow channel (pure fluid), there is no structure available to support other components of the fuel cell, such as the membrane. Methods of the present technology include setting the minimum inverse permeability to a non-zero value and obtaining a grayscale design and fluid velocity field. FIG. 8 provides an exemplary grayscale design representing a porous medium with the inverse permeability κ_(min) being equal to one third of the maximum, for example, κ_(min)=0.333 κ_(max), so there is not any area that is a pure fluid, instead it is all porous. The grayscale of FIG. 8 is now representative of the porous medium fluid flow from the inlet 24 via the tailored permeability fluid streamlines, generally represented by reference number 38, that essentially exits via one outlet.

In various aspects, an effective medium theory approach may be used in order to obtain a width/channel spacing of the microchannels. For example, a functionally graded porous flow medium can be translated into a microchannel network, where the permeability of the porous medium contains embedded information about microchannel dimensions.

Traditional topology optimization does not retain microchannels. However, the present technology utilizes information embedded in the inverse permeability field to guide the relative design of the microchannel structure in order to obtain fine scale features. For example, the inverse permeability can provide a relationship between porous media and microchannel design utilizing an effective medium approach where, for example, the microchannel spacing scales with the square of the channel width. It should be noted that other scale functions may be contemplated. In this regard, the full Navier-Stokes approach is employed for narrow channels for the fluid flow assessments. Each microchannel can be modeled as a fluid-saturated porous medium, numerically solved using the modified Forchheimer-Brinkman-extended Darcy equation for the fluid flow:

${{{{- \frac{d}{dx}}\left\langle p \right\rangle_{f}} + {\mu_{f}\frac{d^{2}}{{dy}^{2}}\left\langle u \right\rangle_{f}} - \frac{\mu_{f}}{K}} \in \left\langle u \right\rangle_{f}} = 0$ which leads to the following relationships:

${\epsilon = \frac{w_{c}}{w}},{K = \frac{\epsilon\; w_{c}^{2}}{12}},{k_{se} = {\left( {1 - \epsilon} \right)k_{s}}},{k_{fe} = {\epsilon\; k_{f}}}$ See, e.g., S. J. Kim, D. Kim, Forced convection in microstructures for electronic equipment cooling, ASME J. Heat Transfer 121 (1999) 635-645. The approach can be utilized to provide a relationship of permeability K as a function of the microchannel width w_(c) and porosity ε, where k_(se), k_(s), k_(fe), and k_(f) are the effective thermal conductivity of solid, thermal conductivity of the solid, effective thermal conductivity of fluid, and thermal conductivity of fluid, respectively. In this regard, the microchannel pitch/spacing is logically linked to the permeability of the design porous material layout; for example, it directly relates the microstructure of the porous material to the dimensions of a channel structure.

Once the grayscale design and fluid velocity fields are obtained, in various aspects, exemplary methods include using the Gray-Scott reaction diffusion equations in order to model and obtain a microchannel layout, for example, to obtain generate a variable pitch microchannel pattern. Once a microchannel layout is obtained, it can be used as a pattern for a flow field area of a bipolar plate.

As is known in the art, reaction-diffusion systems are mathematical models that correspond to physical phenomena. In one example, a change in space and time of the concentration of one or more chemical substances is modeled. In mathematical terms, reaction-diffusion systems generally take the form of semi-linear parabolic partial differential equations represented by the general form as follows: ∂_(t) q=D∇ ² q+R(q), where q(x, t) represents the unknown vector function, D is a diagonal matrix of diffusion coefficients, and R accounts for all local reactions. Reaction and diffusion of chemical species can produce a variety of patterns, reminiscent of those often seen in nature. The Gray-Scott equations model such a reaction and are generally provided as follows:

${\frac{\partial u}{\partial t} = {{D_{u}{\nabla^{2}u}} - {uv}^{2} + {F\left( {1 - u} \right)}}},{\frac{\partial v}{\partial t} = {{D_{v}{\nabla^{2}v}} + {uv}^{2} - {\left( {F + k} \right){v.}}}}$ The partial differential equations modeling this process may be simulated with a variety of numerical techniques. In various aspects, good results can be obtained using methods such as forward Euler integration of the finite-difference equations that one obtains by spatial discretization of the Laplacian, or the diffusion coefficients can be estimated for a target pitch/spacing. For example, in various aspects the value of (pitch){circumflex over ( )}2 (i.e. the square of the microchannel width) can be used to approximate the diffusion coefficient.

The Gray-Scott reaction diffusion model is a mathematical model that describes the behavior of two chemical substances and calculates the concentration of the two substances at a given time based upon the substances diffusion, feed rate, removal rate, and a reaction between the two. This simulation not only models the underlying process of a chemical reaction but can also result in patterns of the substances that are remarkably similar to patterns found in nature. Examples include patterns on animals, such as stripes on zebras, a leopard's skin, spots on butterflies, patterns on fish; fingerprints; ripples on sand; patterns of veins on a leaf; and various other biological phenomena. With the present technology, certain of the patterns resulting from this model can be used for the design of at least a portion of a flow field of a bipolar plate.

To illustrate the model, one analogy presented is to imagine an area or space containing various concentrations of each chemical substance U and V at time zero. Over time, substance U is fed into the reaction at a given rate, while substance V is removed at a given rate. Further, two molecules of V can react with one of U, which converts the substance of U to V as follows: U+2V→3V V→P U, V, and P are chemical substances. A simulation is accomplished using the two partial differential equations listed above, each representing the change in concentration of a substance over time, where u and v are independent variables that represent their respective concentrations; D_(u) and D_(v) are their respective diffusion rates or diffusion tensors, which can be calculated from permeability. The parameter k represents the rate of conversion of V to P; and F represents the rate of process that feeds U and drains U, V, and P. The parameters k and F are arbitrary positive numbers that can be adjusted. Each diffusion tensor is generally a 3×3 matrix reflecting diffusion rates in different directions.

The change in u (upper partial differential equation) is dependent upon its reaction with v (hence the subtraction (−) of uv²) and is fed at a certain rate (+F, scaled to its current concentration). The change in v (lower partial differential equation) is dependent upon its reaction with u (hence the addition (+) of uv²), and is removed at a given rate (−k, scaled by the feed rate and concentration of V). The concentration of U or V at each position is updated at each time increment (typically 1) based upon the result of the corresponding equation. The values for the feed rate, removal rate, and diffusion rate are entered into the equations. For example, on a 2D grid, the Laplacian Operator could be calculated on a convolution matrix. To calculate the new concentration, the current concentration and each surrounding concentration is multiplied by the corresponding value in the matrix (where the current position corresponds to the center position in the convolution matrix) and all values summed. This value technically represents the difference in concentrations between the current position and the surrounding positions.

The above model can be programmed in a suitable computer code as is known in the art. In various aspects, the resulting model provides an image representing the reaction container, with each point or pixel of the image representing the concentration of V (v) at that position. For example, the Gray-Scott reaction diffusion equations can be solved with various numerical methods, such as the finite differential method or the finite element method. The initial value of u and v could be random noise distribution. By solving the Gray-Scott equations for some time period with an appropriate parameter set, a steady state is obtained. Then, the obtained distribution is interpreted as a channel structure. For instance, u>th may represent the wall domain and u<th may represent the channel domain, where th is a threshold value.

Notably, the diffusion tensor is coupled with the permeability distribution of the flow system to be designed. For example, the permeability distribution is projected to the channel pattern with a variable array pitch through the diffusion coefficient distribution of the Gray-Scott equations. The circles A and B in FIG. 8 illustrate, by way of non-limiting examples, certain areas of the flow field that may be provided with a fine pitch, represented by circle A, and with a coarse pitch, represented by circle B.

In various aspects, the present technology may use the Gray-Scott reaction diffusion equations with the grayscale design and fluid velocity field to obtain an anisotropic microchannel layout. This may include using an anisotropic diffusion tensor with the Gray-Scott equations. FIG. 9 illustrates an exemplary method flow chart for designing a pattern for an inlet region of a flow field using an anisotropic diffusion tensor.

In various aspects, a weak and strong anisotropic tensor set is alternatively used to obtain the connected line and space pattern, along with the flow direction. With an isotropic diffusion tensor, the ratio of the microchannel pitch (a) to the length (b) can be represented as a:b=1:1. For a weakly anisotropic tensor, the ratio of a:b is between about 1:1 to about 1:5. For a strongly anisotropic tensor, the ratio of a:b may be up to about 1:100. As noted in FIG. 9, the lateral direction defines the channel pitch. By stretching the tensor in the longitudinal direction, which is the flow direction, the line and space pattern can be obtained. In various aspects, at least two sets of diffusion tensors are used to repeatedly solve the reaction diffusion equations: an isotropic or weakly anisotropic set, and a strongly isotropic set. In various methods, the Gray-Scott equations may be solved for some time period, where the diffusion tensor sets alternate in strength from weak to strong. For example, the diffusion tensor may be switched from isotropic→anisotropic→isotropic. This may be repeated several times, preferably finishing with either a weakly anisotropic set or completely isotropic set. FIG. 10 illustrates an exemplary schematic generally summarizing methods for designing a pattern for an inlet region of a flow field.

Pending patent application Ser. No. 16/275,394, filed on Feb. 14, 2019, which has the same inventors and common ownership with the present application, provides various illustrations showing the differences in design of the flow field based on the type of diffusion tensor(s) used, and is incorporated by reference herein in its entirety. As provided therein, an initial random pattern can be obtained from the Gray-Scott equations using a single diffusion tensor to obtain an inlet flow field using a single isotropic diffusion tensor set, with a ratio of a:b as 1:1. In contrast, an inlet flow field using a single weakly anisotropic diffusion tensor set can be obtained from the Gray-Scott equations, with a ratio of a:b as 1:5, and an inlet flow field using a single strongly anisotropic diffusion tensor set obtained from the Gray-Scott equations, with a ratio of a:b as 1:100. Inlet regions of flow fields may be provided with a combination of diffusion tensors resulting in an optimized microchannel design, including: alternating from strongly anisotropic to isotropic; alternating from isotropic, to strongly anisotropic, to isotropic; alternating from strongly anisotropic, to isotropic, to strongly anisotropic, back to isotropic.

Once the pattern of the flow field is generally designed with a microchannel layout, various channels of the flow field, at least in the inlet region 32, are then defined having discrete areas of discontinuity (reference number 58 of FIGS. 15, 16, and 18) that ultimately provide an increased uniformity in fluid flow. In various aspects, the discrete areas of discontinuity 58 are aligned and located along or near spaced apart contour lines 54 that may be based on a predetermined performance parameter of the bipolar plate. For example, the predetermined performance parameter may be at least one of pressure flow, discussed below with respect to FIGS. 11-16, and velocity flow, discussed below with respect to FIGS. 17-19.

FIG. 11 illustrates an exemplary pressure field 50 obtained from an optimized porous media flow field model for flow from an inlet 24 and through the inlet region 32, exiting an outlet window region 52 ultimately leading through to the other side of the flow field 28 (FIG. 1). As shown, a plurality of spaced apart contour lines 54 are provided substantially aligned with predetermined boundaries of changes in pressure. The contour lines 54 indicate pressure boundaries as illustrated in a color range from dark red (highest pressure) to dark blue (lowest pressure), with the color scale providing contour pressure (in units of Pa) and the numerical scale providing contour pressure (in units of Pa). FIG. 12 illustrates a line graph of a velocity profile associated with the pressure field of the optimized porous media inlet flow field model of FIG. 11 that corresponds to one theoretical outlet flow channel. For purposes of providing an example flow field, FIGS. 13-19 are intended to have 113 representative flow channels at the outlet window region 52 and extending in a direction along the x-axis (FIG. 1).

FIG. 13 illustrates a resulting pressure field 50 obtained from a continuous channel design of the optimized porous media inlet flow field model of FIG. 11. As noted in the area designated with reference number 56, a rate of branching of the channels 40 is not consistent, leading to pressure and flow non-uniformity. As shown, a plurality of spaced apart contour lines 54 are provided substantially aligned with predetermined boundaries of changes in pressure. The contour lines 54 indicate pressure boundaries as illustrated in a color range from dark red (highest pressure) to dark blue (lowest pressure), with the color scale providing contour pressure (in units of Pa) and the numerical scale providing contour pressure (in units of Pa). The pressure field 50 of FIG. 13 indicates a non-uniformity in pressure at the outlet window region 52. FIG. 14 illustrates a line graph of a velocity profile associated with the pressure field of the continuous channel inlet flow field design of FIG. 13, clearly indicating the flow non-uniformity between the flow channels.

FIG. 15 illustrates a resulting pressure field 50 obtained from a discontinuous channel design of the optimized porous media inlet flow field model of FIG. 11. The pressure field 50 is provided having a substantially uniform inlet feed rate in the region 60 immediately adjacent the inlet 24. The pattern includes a plurality of channels 40 defined by a plurality of side walls 42 with discrete areas of discontinuity 58. The contour lines 54 indicate pressure boundaries as illustrated in a color range from dark red (highest pressure) to dark blue (lowest pressure), with the color scale providing contour pressure (in units of Pa) and the numerical scale providing contour pressure (in units of Pa). The pressure field 50 of FIG. 15 having the discrete areas of discontinuity 58 indicates a uniformity in pressure at the outlet window region 52. As shown, the discontinuous channel design of the present technology is able to nearly exactly recover the porous model pressure field distribution provided in FIG. 11. In various aspects, locating the discrete areas of discontinuity 58 with the plurality of contour lines 54 of the pressure field provides a variation of fluid flow pressure between channels of the outlet window region 52 of the inlet region 32 of the flow field and the porous media model domain of less than about 7%, less than about 5%, and even less than about 3%.

FIG. 16 illustrates an exemplary design of a flow field 62 with a discontinuous channel inlet flow field design according to another aspect of the present technology. Differentiated from FIG. 15, the design of FIG. 16 provides a pattern having a non-uniform inlet feed rate in the region 64 immediately adjacent the inlet 24. The pattern includes a plurality of channels 40 defined by a plurality of side walls 42 with discrete areas of discontinuity 58. As shown, many of the discrete areas of discontinuity 58 are aligned and located with contour lines 54 representative of a predetermined performance factor of the bipolar plate. The flow field 62 of FIG. 16 having the discrete areas of discontinuity 58 similarly indicates a uniformity in pressure at the outlet window region 52.

FIG. 17 illustrates an exemplary velocity field 66 obtained from an optimized porous media inlet flow field model for flow from an inlet 24 and through the inlet region 32, exiting an outlet window region 52 ultimately leading through to the other side of the flow field 28. As shown, a plurality of spaced apart arrows 68 are provided to represent direction and relative magnitude of the velocity. The contoured color regions indicate a central area 70 with a higher velocity and indicate velocity boundaries as illustrated in a color range from dark red (highest velocity) to dark blue (lowest velocity), with the color scale providing contour velocity (in units of m/s).

FIG. 18 illustrates a resulting velocity field 66 obtained from a discontinuous channel design model of the optimized porous media inlet flow field model of FIG. 17. The pattern includes a plurality of channels 40 defined by a plurality of side walls 42 with discrete areas of discontinuity 58. The contoured color regions indicate velocity boundaries as illustrated in a color range from dark red (highest velocity) to dark blue (lowest velocity), with the color scale providing flow velocity (in units of m/s). The velocity field 66 of FIG. 18 having the discrete areas of discontinuity 58 indicates a uniformity in flow velocity at the outlet window region 52. As shown, the discontinuous channel design of the present technology is able to nearly exactly recover the porous model velocity field distribution provided in FIG. 17. In various aspects, locating the discrete areas of discontinuity 58 with the plurality of contour lines of the velocity field provides a variation of fluid flow velocity between channels of the outlet window region 52 of the inlet region 32 of the flow field and the porous media model domain of less than about 7%, less than about 5%, and even less than about 3%.

FIG. 19 is a graph that shows the results of a comparison between a discontinuous channel inlet velocity field and a target specification velocity field. The average outlet velocity is 0.124 m/s. The average flow non-uniformity is about 4.5%, relative to the average outlet velocity. The maximum flow non-uniformity is about 9.7%, relative to the average outlet velocity.

Once an optimized microchannel layout and design is obtained and incorporated into a pattern, the present technology also provides for the manufacturing of the channels and creating the flow field of the bipolar plate, optionally with a tailored porosity. FIG. 20 provides a cross-section view of an exemplary set of basic channels 40, walls 42, and bottom wall 44 of an inlet region 32 of a flow field that can be made according to various manufacturing techniques. It should be understood that the specific design of the channels and walls can be quite detailed in various aspects.

Numerous materials and fabrication techniques can be used, in combination with the pattern, to create a network of a plurality of channels/microchannels 40 defined by a plurality of side walls 42 in various regions of a flow field. In various aspects, the present technology is particularly useful with functionally graded materials and/or microchannels, a porous media, or a hybrid combination thereof. A functionally graded material (“FGM”), also referred to as a functionally graded porous material, is generally known in the art as a material with a changing composition, microstructure, and/or porosity across a volume of the material, in one or more direction. FGMs can be specifically designed to perform a set of specified functions, and have almost endless possibilities for tailoring and customization, which can be very useful in making a flow field, or portion thereof, for a bipolar plate. Design parameters of the FGM can be used with the effective medium theory, as described above, in order to translate a FGM into a microchannel network, where at least one dimension of the microchannels, such as a width, is based on a permeability of the FGM. FGMs can generally be classified as being porosity and pore size gradient-structured, chemical/composition gradient-structured, and microstructural-gradient structured; each of which can be used singly or in combination with the flow field patterns designed herein. Not only does the present technology envision shaping and/or machining an existing FGM into a pattern, but also contemplates the creation of FGMs using additive manufacturing and similar techniques and technologies. For example, microstructural gradient functionally gradient material can be additive manufactured such that the microstructure is tailored so that different microstructures are produced in the material, which is made to change gradually, so as to achieve the desired properties from the material. In certain aspects, microstructural gradation can also be achieved during a solidification process using quenching. In other aspects, sintering techniques can be useful to produce a varying microstructure within a material.

In the broadest sense, additive manufacturing techniques can be used to make a pattern in the flow field, or portion thereof, where it is desired that at least a portion of the flow field is custom made and tailored to a specific design. Various non-limiting examples of additive manufacturing techniques include fabricating using laser sintering techniques, selective laser melting and electron-beam melting techniques, 3-D printing, and/or using binder jet technology that may use a combination of printing and curing. In other aspects, various electroplating techniques may also be used for fabricating or working with a functionally graded porous medium. Certain composite materials, polymers, thermoplastics, resins, and the like may be considered as useful materials, depending on the patterns and designs.

In still other aspects, metal inverse opal (“MIO”) technology can be used with a suitable metal, such as copper, nickel, titanium, and the like. MIO technology can be tailored and customized as is known in the art. In one common example, a functionally graded porous medium can generally be obtained by templating polystyrene particles of different sizes, and/or controlling a sintering time for different templated regions. If desirable, MIO techniques can be used that create necks, pores, and fluidic routing structures exerting specific capillary forces, wicking, and flow patterns. Conventional metal machining and/or stamping techniques can be used in combination with the above techniques to ultimately provide the bipolar plates having flow fields with variable pitch microchannel designs.

Various hybrid porous structured materials and combinations of materials and methods of manufacture can also be used with the present technology. Hybrid porous structured materials may include a porous region and a non-porous region. Metal foams, and the like, may also be incorporated with designs of the present technology as a functionally graded porous material. In various aspects, dense, non-porous channel walls 42 can be provided with a coating of a porous material. Certain designs may include dense walls 42 with a porous bottom wall 44. In other aspects, the channel walls 42 may be made of a porous material or FGM, and the bottom wall 44 can be non-porous. The changes in porosity can also be gradient. In still other aspects, one or more regions of the bottom wall 44 can be designed to have a first degree of porosity, while the channel walls 42 may have a second degree of porosity, and so on. In certain aspects, the bottom wall 44 may be a part of the substrate or bipolar plate, and a plurality of microchannels may be subsequently created from a FGM, or formed on the bipolar plate using additive manufacturing, electroplating, or a combination thereof in order to create a variable pitch microchannel pattern configured to direct fluid from an inlet to an outlet.

The foregoing description is provided for purposes of illustration and description and is in no way intended to limit the disclosure, its application, or uses. It is not intended to be exhaustive or to limit the disclosure. Individual elements or features of a particular embodiment are generally not limited to that particular embodiment, but, where applicable, are interchangeable and can be used in a selected embodiment, even if not specifically shown or described. The same may also be varied in many ways. Such variations should not be regarded as a departure from the disclosure, and all such modifications are intended to be included within the scope of the disclosure.

As used herein, the phrase at least one of A, B, and C should be construed to mean a logical (A or B or C), using a non-exclusive logical “or.” It should be understood that the various steps within a method may be executed in different order without altering the principles of the present disclosure. Disclosure of ranges includes disclosure of all ranges and subdivided ranges within the entire range, including the endpoints.

The headings (such as “Background” and “Summary”) and sub-headings used herein are intended only for general organization of topics within the present disclosure, and are not intended to limit the disclosure of the technology or any aspect thereof. The recitation of multiple embodiments having stated features is not intended to exclude other embodiments having additional features, or other embodiments incorporating different combinations of the stated features.

As used herein, the terms “comprise” and “include” and their variants are intended to be non-limiting, such that recitation of items in succession or a list is not to the exclusion of other like items that may also be useful in the devices and methods of this technology. Similarly, the terms “can” and “may” and their variants are intended to be non-limiting, such that recitation that an embodiment can or may comprise certain elements or features does not exclude other embodiments of the present technology that do not contain those elements or features.

The broad teachings of the present disclosure can be implemented in a variety of forms. Therefore, while this disclosure includes particular examples, the true scope of the disclosure should not be so limited since other modifications will become apparent to the skilled practitioner upon a study of the specification and the following claims. Reference herein to one aspect, or various aspects means that a particular feature, structure, or characteristic described in connection with an embodiment or particular system is included in at least one embodiment or aspect. The appearances of the phrase “in one aspect” (or variations thereof) are not necessarily referring to the same aspect or embodiment. It should be also understood that the various method steps discussed herein do not have to be carried out in the same order as depicted, and not each method step is required in each aspect or embodiment. 

What is claimed is:
 1. A method for designing a microchannel layout for a flow field of a bipolar plate, the method comprising: defining a fluid flow optimization domain with boundary conditions and loads, the fluid flow optimization domain including an x-axis defining a longitudinal direction from an inlet leading to an outlet of the flow field, a y-axis defining a transverse direction with respect to the x-axis, and a z-axis perpendicular to both the x-axis and the y-axis; using a gradient-based algorithm together with computational fluid dynamics to optimize a porous media model domain; setting a minimum inverse permeability to a non-zero value, and obtaining a grayscale design and fluid velocity field; using Gray-Scott reaction diffusion equations with the grayscale design and fluid velocity field to obtain a microchannel layout with a plurality of channels; providing the plurality of channels with discrete areas of discontinuity; and incorporating the microchannel layout as a pattern for an inlet region of the flow field of the bipolar plate.
 2. The method according to claim 1, further comprising developing a plurality of contour lines in the flow field, wherein the plurality of contour lines are spaced apart based on a predetermined performance parameter of the bipolar plate.
 3. The method according to claim 2, further comprising locating the discrete areas of discontinuity with the plurality of contour lines.
 4. The method according to claim 3, further comprising determining an optimized pressure field for the flow field, wherein the optimized pressure field is used as the predetermined performance parameter to define the plurality of contour lines based on differences in pressure.
 5. The method according to claim 4, wherein locating the discrete areas of discontinuity with the plurality of contour lines of the pressure field provides a variation of fluid flow pressure between channels of the outlet of the flow field and the porous media model domain of less than about 5%.
 6. The method according to claim 3, further comprising determining an optimized velocity field for the flow field, wherein the optimized velocity field is used as the predetermined performance parameter to define the plurality of contour lines based on differences in velocity.
 7. The method according to claim 6, wherein locating the discrete areas of discontinuity with the contour lines of the optimized velocity field provides a variation of fluid flow velocity between channels of the outlet of the flow field and the porous media model domain of less than about 5%.
 8. The method according to claim 1, further comprising using the gradient-based algorithm together with computational fluid dynamics to optimize the porous media model domain for minimum flow resistance across the porous media model domain.
 9. The method according to claim 8, wherein using the gradient-based algorithm together with computational fluid dynamics to optimize the porous media model domain for minimum flow resistance comprises using a modified Darcy flow equation for porous media.
 10. The method according to claim 1, further comprising using the gradient-based algorithm together with computational fluid dynamics to optimize the porous media model domain for a uniform flow velocity in a direction of the x-axis.
 11. The method according to claim 1, further comprising using the gradient-based algorithm together with computational fluid dynamics to optimize the porous media model domain for a minimum flow velocity in a direction of the y-axis.
 12. The method according to claim 1, further comprising using Gray-Scott reaction diffusion equations with the grayscale design and fluid velocity field to obtain an anisotropic microchannel layout with a plurality of channels.
 13. The method according to claim 12, further comprising using a weak and strong anisotropic diffusion tensor set alternatively to obtain a connected line and space pattern for the anisotropic microchannel layout.
 14. A method for designing an anisotropic microchannel layout for a flow field of a bipolar plate, the method comprising: defining a fluid flow optimization domain with boundary conditions and loads; using a gradient-based algorithm together with computational fluid dynamics to optimize a porous media model domain; setting a minimum inverse permeability to a non-zero value, and obtaining a grayscale design and fluid velocity field; using Gray-Scott reaction diffusion equations and an anisotropic diffusion tensor along with the grayscale design and fluid velocity field and to generate an anisotropic microchannel layout having a connected line and space pattern including a plurality of channels; providing the plurality of channels with discrete areas of discontinuity; and incorporating the anisotropic microchannel layout as a pattern for an inlet region of the flow field of the bipolar plate to provide flow uniformity.
 15. The method according to claim 14, further comprising using the gradient-based algorithm together with computational fluid dynamics to optimize the porous media model domain for a minimum flow resistance across the domain.
 16. The method according to claim 14, further comprising using the gradient-based algorithm together with computational fluid dynamics to optimize the porous media model domain for a uniform velocity across a direction of the domain.
 17. The method according to claim 14, further comprising using the gradient-based algorithm together with computational fluid dynamics to optimize the porous media model domain for a minimum velocity across a direction of the domain.
 18. The method according to claim 14, further comprising using a weak and strong anisotropic diffusion tensor set alternatively to solve the Gray-Scott reaction diffusion equations, and determining a spacing of microchannels based on effective medium theory, wherein the discrete areas of discontinuity are disposed along spaced apart contour lines based on a predetermined performance parameter of the bipolar plate selected from at least one of pressure flow and velocity flow.
 19. A bipolar plate for a fuel cell, the bipolar plate comprising: an inlet; an outlet; and a flow field having an inlet region comprising a pattern defining a plurality of anisotropic microchannels configured to provide fluid communication between the inlet and the outlet, the anisotropic microchannels defined having discrete areas of discontinuity, wherein the pattern is designed using an inverse permeability field and solving Gray-Scott reaction-diffusion equations with an anisotropic diffusion tensor to obtain an anisotropic microchannel layout with channel spacing based on effective medium theory.
 20. The bipolar plate according to claim 19, wherein the discrete areas of discontinuity are disposed along spaced apart contour lines based on a predetermined performance parameter of the bipolar plate selected from at least one of pressure flow and velocity flow. 